Measuring the condensate fraction of rapidly rotating trapped boson systems: 
off-diagonal order from the density profile 
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We demonstrate a direct connection between the density profile of a system of ultra-cold trapped 
bosonic particles in the rapid-rotation limit and its condensate fraction. This connection can be 
used to probe the crossover from condensed vortex-lattice states to uncondensed quantum fluid 
states that occurs in rapidly rotating boson systems as the particle density decreases or the rotation 
frequency increases. We illustrate our proposal with a series of examples, including ones based on 
models of realistic finite trap systems, and comment on its application to freely expanding boson 
density profile measurements. 
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When a system of coherent trapped bosons is brought 
to equihbrium in a rotating frame of reference, vortices 
are induced in the order parameter of its condensate and 
centrifugal forces lead to radial expansion. The influ- 
ence of the rotating frame may be described mathemat- 
ically by introducing an effective magnetic field oriented 
along the axis of rotation and reducing the radial confine- 
ment strength. Taking this point of view, very rapidly 
rotating boson systems enter the quasi-two-dimensional 
(2D) strong-field quantum Hall regime in which inter- 
actions are always important. Many-fermion systems in 
the quantum Hall regime have been thoroughly studied 
over the past twenty years and have exhibited a wide va- 
riety of subtle and novel behaviors, many of which were 
completely unanticipated by theory . This history pro- 
vides a powerful motivation for creating boson quantum- 
Hall systems. Most current experiments in rotating Bose- 
Einstein systems create condensates in which the effec- 
tive magnetic field is relatively weak and axial-direction 
quantum confinement effects are negligible when com- 
pared to interaction strengths. For such condensates, 
the nucleation of vortices and the formation of vortex 
lattices, together with their decay and evolution, have 
been observed experimentally and studied theoretically 
by several groups ^^^^^Jt,^. 

The limit of quasi-2D bosons in the rapid-rotation (or 
equivalently, quantum Hall) regime has so far been stud- 
ied only theoretically 0, [3, E E El • The physics 
of trapped bosons in the quantum Hall regime is enriched 
by the near-degeneracy of single-particle energy levels, 
which enhances the importance of condensate quantum 
fiuctuations 0, reduces the condensate fraction, and 
eventually at low boson densities, leads to quantum fluid 
states nahj that are expected to have unusual proper- 
ties analogous to those that have been uncovered in stud- 
ies of the fermion quantum Hall effect. One of the major 
stumbling blocks in contemplating studies of bosons in 
the quantum Hall regime has been the lack of a reliable 
tool to directly measure the degree of Bose- Einstein con- 
densation and the accompanying long-range many-boson 



phase coherence. Unlike its familiar zero-rotation cousin, 
a vortex-lattice condensate in the quantum Hall regime 
does not have a characteristic velocity distribution that 
allows it to be identified from the velocity distribution 
of the released system. In this paper, we show that in 
the quasi-2D rapid-rotation limit there is, however, a sur- 
prisingly simple relationship between diagonal elements 
of the many-boson density-matrix and the off-diagonal 
order responsible for coherence. The condensate fraction 
can be probed simply by measuring the spatial distribu- 
tion of the boson density. 

Before explaining the simple but powerful relation that 
we propose be exploited, we first discuss its regime of 
validity. A system of N interacting bosons of mass M in 
a cylindrical trap (with radial and axial trap frequencies 
rir and Qz) that is rotating with angular velocity Qz is 
well-described by the rotating-frame Hamiltonian 15] 
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where V{r)=gS{r) is a hard-core interaction potential of 
strength g — Anh^as/M , and is the s-wave scatter- 
ing length. This Hamiltonian is equivalent to that of 
a system of charge- Q bosons with weakened radial con- 
finement, under the influence of a magnetic fleld B = 
{2Mn/Q)z, with cyclotron frequency fi^ = QB/M = 2^1 
and magnetic length i — yjhj (2A/il). The Hilbert space 
of this system is spanned by the axial (z-direction) quan- 
tum states and the various Landau-level states produced 
by the effective magnetic field. For close to $7,-, the ef- 
fective confinement potential becomes very weak, which 
causes the bosons to spread out radially and, in the limit 
of interest to us here (see below), allows all bosons to be 
accomodated in the lowest-energy axial quantum state. 

In the quasi-2D limit, the boson system can be 
described by a 2D model whose short-range interac- 



tions have strength g2D 
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^JhjMVTz- The areal density of bosons in such a trap 
can be estimated using a 2D Thomas-Fermi (TF) ap- 
proximation in which Mffi? — r2^)r^/2 + g^Dn-iD is con- 
stant inside the system fl(j| Since the bosons can form 
rather uniform-density clouds (of the type favored by 
short-range interactions) entirely within the LLL, there is 
little incentive for them to occupy higher kinetic-energy 
Landau levels, unless the typical interaction strength 
920^20 is much larger than the Landau-level splitting 
2MI. The TF approximation boson Landau- level fill- 
ing factor is given by v(r) = 2'K(.'^n2D{f) — — 
r^lB?), where = ^20^ It^M{VlI - 1^2), and vl = 
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characterizes the degree of Landau-level mixing induced 
by interactions. Neglecting boson Landau-level mixing, 
as we do in the rest of the paper, is a good approxima- 
tion provided that is not large. Similarly, neglecting 
mixing of higher axial energy eigenstates is a good ap- 
proximation provided that 7^ = 2jLLn/nz is not too 
large. In mean-field theory, the effect of a finite value 
of jz is to weaken axial confinement at the center of the 
trap and therefore to effectively weaken g2D- The conclu- 
sions we reach below would not be altered by a position- 
dependent effective interaction strength. We therefore 
expect our results to apply for values oi^LL and 7z that 
are substantially larger than one, although more detailed 
calculations will be required to establish the range of the 
rapid-rotation limit more quantitatively. Values of 7x,l 
and 72 that are not much larger than unity are already 
within reach of current experimental designs [l7j. 

Condensation of a system of trapped bosons (we con- 
sider here only the zero-temperature case) is discussed 
most generally in terms of the spectrum of its one-particle 
density matrix (ODM), whose position representation is 
related to the many-boson wave function by 



pi(r,r') = 7V / dr2---dYN^l,{T,T2, 
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x*Ar(r',r2,-- - ,rAr) = (^t(r)^(r')), (3) 

where pi(r,r) = n(r) is the boson density, and the trace 
of the ODM is the total particle number. A''. If one or 
more of the eigenvalues, Ai, of the ODM is extensive (i.e., 
Xi oc N), then the system is Bose-condensed In a 

particle-conserving system, the ratio 
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will be of order unity only if the system is Bose-condensed 
and the condensate fraction will be y^, up to a correction 
of order 

The many-body physics of quasi-2D boson systems in 
a magnetic field simplifies in the rapid-rotation limit be- 
cause the many-boson Hilbert space can be projected 



to the lowest Landau level (LLL). These simplifications 
are often most conveniently captured in the symmet- 
ric gauge, where they imply that many-particle wave 
functions must have the form \E'7v(ri, r2, • • • ,ym) = 
J[z\ exp (— |zfe|2/4£2) ^ where f[z\ is analytic in each 
of the N complex bosonic coordinates = Xk +iyk- The 
analyticity property has played a key role in achieving 
an understanding of many-fermion physics in the quan- 
tum Hall regime, for example in Laughlin's recognition 
is] that incompressible states would occur at certain 
Landau-level filling factors. Because of this property, the 
ODM is completely specified by the functional form of 
its diagonal matrix element. For fermions, this property 
leads [13 to Hartree-Fock exchange energies that depend 
only on the particle density; for bosons, it has more fun- 
damental implications, as we now discuss. 
It follows from of \I'jv(ri,r2, • • • ^vn) that 

pi(r', r) = n[z\z'] exp (^-^^^ - + ^ j , (5) 

where n[z* ^ z] is the boson density expressed as a function 
of z = X + iy and z* — x — iy. Equation js]) can be used 
to express (3 in terms of the 2D Fourier transform of the 
boson density: ?i[r] — n[z*,z] — J dqn{q) cxp[i{q* z + 
gz*)/2]/(47r2) , where and q = Qx + iqy We find that 
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We propose that Eq. lO be used to measure the conden- 
sate fraction of rapidly rotating quasi-2D Bose-Einstein 
condensates. 

To understand the strengths and weaknesses of this 
condensate measure as an experimental probe, it is help- 
ful to discuss several different examples. We first consider 
the case of a rapidly rotating fully condensed BEG that 
contains a single vortex; for typical traps and interaction 
strengths, such a state would occur only for boson par- 
ticle numbers that are too small to study at the present 
time. In this case, the Gross-Pitaevskii (GP) order pa- 
rameter m is (j>{z) = .yN/{4:n£^)zexp{-\z\^/4f) which 
leads to n{q) = N{1 - q'^i'^ /2) cxp{-q^f /2), and yields 
/3 = 1, as expected. 

A qualitatively similar but more realistic case is illus- 
trated in Fig.n The upper panel of this figure shows the 
mean-field-theory 2D boson density profile, obtained by 
solving the GP equations following Refs. 0and|23, while 
allowing the trap's rotational symmetry to be broken and 
optimizing with respect to the number of vortices, for 
the parameters flr/fl = 1.04, g2D/(?i^M) = 0.0058, and 
N — 10, 000. These parameters were chosen to closely 
match those of current experimental systems ^21]. The 
lower graph in Fig. ^ shows the profile of the angle- 
averaged boson filling factor. Under these circumstances, 
we expect that the condensate fraction will be large, 
at least in the limit of zero temperature. When smoothed 
radially over a length ^ £, the boson density profile 
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FIG. 1: Upper figure: Density profile for the ground state of a 
rapidly rotating boson system with Qr/^ ~ 1.04, g2DM/f? — 
0.0058, and A'^ = 10, 000. These parameters are similar to 
ones that can be achieved by current experimental systems. 
For these parameters, the optimal number of vortices is 19. 
Lower figure: Angle-averaged boson filling-factor profile with 
radial broadening Ar. The long-dashed line in this figure is 
the TF approximation boson filling-factor profile for the same 
parameters. 



agrees reasonably well with that predicted by TF the- 
ory, illustrating the dominant role of interactions in de- 
termining the course-grained density profile, even in the 
rapid-rotation limit. In Fig. El we show the integrand 
of Eq. © that is involved in calculating /3, for the con- 
figuration shown in Fig. ^ In practice however, only 
the first few rings of Fig. [21 can be resolved in current 
experiments and a finite wave-vector cutoff must be in- 
troduced |13]. In this case, two different estimates of 
(3 might prove useful. A lower bound /3_ is obviously 
obtained by calculating (5 using an experimentally con- 
strained momentum cutoff Qmax- A separate estimate of 
/3, which is more likely to be useful in practice for current 
experimental systems, is motivated by observing that 



FIG. 2: |n(g)|^ exp(g^^^/2) for the density profile shown in 
Fig. The axes are those of the 2D wave vector, multiplied 
by the effective magnetic length, i. 



in the case of a quasi-2D (untrapped) condensate each 
Brillouin-zone in momentum space contributes equally 
to / d^q|n(q)p exp(|(j2|^2^2). It follows that the ratio of 
diffraction peak strengths integrated over Brillouin-zones 
centered on q (condensate contributions only) and 
q = (condensate plus normal fluid contributions) can 
be used to estimate (3. The efficacy of this approach 
could be checked by comparing estimates from differ- 
ent shells . For the example of Fig. [3 we obtain 
13- w 0.27, 0.57, 0.83 for cutoffs of q^^J = 4, 6, 9, and, 
using the first shell for the second estimate we obtain 
(3 « 0.93; i.e., all show clear evidence of Bose-Einstein 
condensation, albeit with an inaccurate value of the con- 
densate fraction. (The condensate fraction is, of course, 
always one for GP approximation density profiles.) As 
this calculation illustrates, we do not anticipate that our 
method will be especially useful for measuring small de- 
viations of the condensate fraction from unity in strongly 
condensed systems. We do, however, expect that mean- 
ingful measurements of (3 will be possible when the con- 
densate fraction is strongly reduced due to either thermal 
or quantum fluctuations, and anticipate that strong re- 
ductions will occur at available temperatures even at very 
high values of vq. 

We next examine the quantum fluid regime in which 
rotational symmetry is not broken and the 2D density 
profile should be given approximately by the TF approx- 
imation expression, n{r) — vq{1 — / R^) / {2'kP). For 
R/i > 1, we find that (3 « (2/3)i/o/iV, after introduc- 
ing an ultraviolet cutoff to remove an artificial cusp in 
the TF density. This result illustrates the ability of our 
condensate measure to capture the property that, in the 
rapid-rotation limit, BEC can occur only when accom- 
panied by broken rotational symmetry. Incompressible 
boson states ^2 Q| at very small filling factors are 
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expected to produce more uniform density profiles, but 
they also have f3 ~ ly/N and are thus uncondensed. 

As a final example, we consider a state with two in- 
dependent condensates, A and B, with a density profile 
n(r) — XAUAir) + Xsnsi'T), where \a + \b — 1- For the 
case in which species A is in an uncondensed constant- 
density state and species B is a large vortex-lattice con- 
densed state, we obtain (3 = \^^ + {2Xa — / N . There- 
fore, (3 will be of order unity even if only one of the species 
is condensed. 

Another experimental issue that might sometimes arise 
in the measurement of (3 is the determination of the 
boson density profile from free-expansion images. As 
shown in Refs. ^3 and |23j the evolution of the order 
parameter under 2D free expansion (i.e., when the con- 
finement in the z-direction is maintained) is given by 
<^iab(r,t) - eM-i{nL, - ^J)T{t)/h]^,,t{Y/Kt))/\{t) , 
for a system obeying the time-dependent GP equation. 
Here (/)iab is the order parameter viewed in the labora- 
tory frame, and (/>rot(r) is the order parameter before 
the trap has been turned off (obtained from the time- 
independent GP equation in the initial rotating frame), 
\{t) = y/T+T^, and rirTit) = arctan(f]rt). The only 
assumption used to obtain this result, aside from ignor- 
ing interactions with the thermal cloud (which is justified 
by the short expansion times involved), is the validity of 
the time-dependent GP equation. In the quantum Hall 
regime bosons can enter a regime in which quantum fluc- 
tuations render these assumptions invalid, particularly 
when regions of the condensate near the edge that have 
low local filling fractions become important. It is possible 



that such low-density regions could evolve under free ex- 
pansion differently than the higher filling fraction regions 
near the center. Hence, whenever P is no longer of or- 
der unity, the analysis of the boson density profile using 
free-expansion measurements is suspect and one must re- 
sort to non-destructive [i.e. unexpanded) measurements 
of the density profile. Given the large radial expansion 
expected in this regime of interest here, however, it may 
be possible to resolve a sufficiently large number of shells 
of reciprocal space fT.(q) diffraction peaks to apply the 
method we propose. 

In summary, we propose a new experimental tool to 
explore the transition between a vortex-lattice Bose- 
Einstein condensate and an uncondensed quantum-fluid 
state in rapidly rotating boson systems. This new tool 
derives from a deep connection between the density pro- 
file of a system of ultra-cold bosonic particles in the rapid- 
rotation limit and their degree of condensation. It applies 
in the quasi-2D rapid-rotation limit that current experi- 
ments are approaching. 
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